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We present a stochastic volatility market model where volatility is correlated with return 
and is represented by an Ornstcin-Uhlcnbcck process. With this model we exactly mea- 
sure the leverage effect and other stylized facts, such as mean reversion, leptokurtosis 
and negative skewness. We also obtain a close analytical expression for the characteristic 
function and study the heavy tails of the probability distribution. 

1. Introduction 

During decades the diffusion process known as the geometric Brownian motion has 
been widely accepted as one of the most universal models for speculative markets. 
It was proposed by Osborne in 1959 1 when he observed that empirical distribu- 
tions of prices were biased and in disagreement with the theoretical distribution of 
Bachelier's arithmetic Brownian dynamics. 

However, specially after the 1987 crash, the geometric Brownian motion and its 
subsequent Black-Scholes (B-S) formula were unable to reproduce the option price 
data of real markets. Several studies have collected empirical option prices in order 
to derive their implied volatility. These tests conclude that the implied volatility 
is not constant and it is well-fitted to a U-shaped function of moneyness whose 
minimum is at moneyness near to 1, i.e., when current stock price is equal to the 
striking price. This effect is known as the smile effect and shows the inadequacy of 
the Black-Scholes model since this assumes a constant volatility 2 . 

A possible way out to this inconsistency is assuming that volatility is not a 
constant but an unknown deterministic function of the underlying price. The de- 
terministic volatility still allows to manage the option pricing within the B-S theory 
although in most of the cases it is not possible to derive an analytic option price. 
Within this approach, there exists the ARCH-GARCH models and their subse- 
quent extensions 3,4 . These models do well in describing the implied volatility but 
their disadvantage is that some of their parameters substantially change with time 
frequency 5 . 

The stochastic volatility (SV) models are another possible choice. These models 
assume the original log-Brownian motion model but, as their name indicates, with 
random volatility. At late eighties several different SV models were presented 6 ' 7 > 8 . 
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All of them propose a two-dimensional process involving two independent variables: 
the stock and the volatility. These works are basically interested on option pricing 
theory and ignore the statistical properties of the market model, although they are 
indeed able to reproduce the smile effect. 

After then there have appeared several papers extending and refining the original 
SV models but again many of them are solely interested on adequately describe 
empirical option prices. Stein and Stein 9 is an exception to this tendency because 
they study the most important statistical properties of a volatility model following 
an Ornstein-Uhlenbeck (O-U) process. 

We believe that the relative small number of works dealing with market dynamics 
based on SV models is essentially due to two reasons: (i) Their statistical properties 
are difficult to be analytically derived, and the analysis is even much more involved 
when there are correlations between volatility and stock, (ii) It is commonly asserted 
that empirical data available are not enough for obtaining a reliable estimation of 
all parameters involved in an SV model 10 . 

Our present work wants to modify these statements for the correlated Ornstein- 
Uhlenbeck stochastic volatility (O-U SV) model because we are able to analytically 
derive the main statistical properties of it. On the other hand, the leverage corre- 
lation recently observed by Bouchaud etal. 11 allows us to estimate all parameters 
involved not only in our O-U SV process but eventually in any SV process. 

Recent research on empirical markets data has provided a set of requirements 
that a good market model must obey. In this regard Engle and Patton 5 have listed 
a number of stylized facts about the volatility. The results we will herein derive 
conclude that the SV models are good candidates fairly accomplishing these stylized 
facts. We will prove this and confront the statistical properties of the correlated 
O-U SV model with the statistical properties of the daily return changes for the 
Dow- Jones stock index. 

The paper is divided in 7 sections. After this introduction, we present our 
stochastic volatility market model in Section 2 and study the statistical properties 
of the volatility in Section 3. Section 4 is specifically focussed on the leverage effect. 
Section 5 is devoted to show how to estimate from the Dow- Jones index (1900-2000) 
the parameters of the model. Finally, Section 6 concentrates on the derivation of the 
probability distribution through obtaining the characteristic function. Conclusions 
are drawn in Section 7 and technical details are left to Appendices. 

2. The stochastic volatility market model 

The starting point of any stochastic volatility model is the log-Brownian stochastic 
differential equation: 

^l=^dt + adW 1 (t), (2.1) 

where /i is the drift and a is the volatility. The SV models refine this dynamics 
taking a = a(t) to be stochastic. There exist a large class of such models but, 
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to our knowledge, the dynamics of a is not definitively associated to any specific 
process (see the monograph of Fouque etal. 10 or the review by Ghysels etal. 12 ). 
We choose the Ornstein-Uhlenbeck (O-U) stochastic volatility model because, as 
we will see shortly, it is one of the simplest approaches still reproducing the main 
observed features of markets. We thus assume that the random dynamics of a(t) is 
given by 9 

da{t) = -a(a - 6)dt + kdW 2 (t). (2.2) 

Equations (2.1) and (2.2) contain a two-dimensional Wiener process (Wi(t), W2(t)), 
where dWi(t) = £,i(t)dt (i = 1,2), and £j(i) is Gaussian white noise processes with 
zero mean, i.e., 

E[£i(t)] = and E[Z i (t)Z j (t')}=p ij S(t-t'). (2.3) 

We note that the cross-correlation is given in terms of the Dirac delta function, i.e., 
S(x) = for all x ^ and 

/ a a \ sv \ i f Mx) — oo < x < a, 
.„ ^ Z)6{X ~ Z)d2 = { otherwise; ^ 

where (f>(x) is an arbitrary integrable function. Note that pij — pji and pa = 1, 
hence the components of p^ are reduced to 

9 ^ A (2 - 5) 

where the parameter p is given by Eq. (2.3), i.e., E [£i(i)£2(i')] = pS(t — t') which, 
in terms of the Wiener process, is equivalent to say that E \dW\{t)dW2{t')] = 
when t ^ t', and 

E[dW 1 {t)dW 2 {t)] = pdt. 

Thus p is the correlation coefficient and it has no definite sign (—1 < p < 1). 
However, it is known that a negative p is able to provide the skewness observed 
in financial markets 10 . One of our objectives here is not only to show that p is 
negative but also to estimate its value from empirical data. 

In what follows it turns to be more convenient to work with the zero-mean return 
defined as 

dX=^- fidt, (2.6) 

and whose SDE reads 

dX(t) = a(t)dW 1 (t). (2.7) 

The zero-mean return X{t) has a simpler dynamics than the stock price S{t) because 
it only contains the random fluctuation <jdW\ . Nevertheless, this process still retains 
the most interesting features of the whole dynamics. In the Appendix A, we give 
en explicit expression for X(t) and derive some key features. 
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3. The Ornstein-Uhlenbeck volatility process 

We will now present the main properties of the 0-U volatility. The starting point 
is the solution of Eq. (2.2): 

a(t) = <Toe- a( *-* o) + 6(1 - e- a{t ~ to) ) + k f e^*"*' 3 >dW 2 (t'), (3.8) 

where we have assumed that the process started at time t = to, when volatility 
was (To. From now on we will assume that the volatility is in the stationary regime. 
This means that the market started long time ago thus t — > — oo and the stationary 
volatility reads 

a(t)=6 + k e- a ^ t - t 'UW 2 (t'), (3.9) 

J — oo 

whose average value, variance and correlation are 

E [a] = 6, Var[cr] = E [a 2 ] - E [a] 2 = k 2 /2a, (3.10) 

and 

E [a(t + r)a(t)} =6 2 + [k 2 /2a)e' aT . (3.11) 

These expressions provide a physical interpretation of the parameters of the model, 
especially 9, the expected volatility, and a, the inverse of the volatility "correlation 
time" (see below). 

Note that the stationary volatility (3.9) is a linear functional of the Wiener 
process. Therefore, it is a Gaussian process (the same applies to "transient" volatil- 
ity (3.8)). In consequence the stationary distribution is uniquely determined by the 
mean and variance given in Eq. (3.10). That is, 

P*fr) = - / ;2/ e-^- B)lk \ (3.12) 
2^Jixk l J a 

Before proceeding further we want to address the question of the sign of <r(t). 
If one adopts the O-U process (2.2) as a model for stochastic volatility one may 
argue that <r(t) has no definite sign which can be seen as an inconvenience for a 
"good" SV model. Let us see that this is not really the case. First of all, the actual 
evaluation of volatility is very difficult, not to say impossible, since volatility itself 
is not observed. In practice, the so-called instantaneous volatility is derived from 

lhnJ[X(t + At)-X(t)] 2 /At, (3.13) 

where we have used the zero-mean return defined above. Due to the limit At — ► 0, 
this equation has to be taken as an infinitesimal difference. Wc thus define 



instantaneous volatility = \J dX{t) 2 / dt. 



(3.14) 
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From Eq. (2.7) and the fact dW 2 = dt, we get 



instantaneous volatility = \J u{t) 2 = \a(t)\ . (3.15) 
Observe that in this definition no sign is attached to the random variable u(t). 
3.1. The correlated process 

As we have mentioned in Section 2, the model permits correlations between the stock 
and the volatility. We will now examine the important effects of these correlations. 
From Eq. (3.9), we see that the correlation between the stationary volatility and 
the random component of return variations, dW\(t), is 

[t + T 

E [a(t + r)dWi(t)] = k e -«{t+r-t') E [dW 2 (t')dW 1 (t)\ , 

J — oo 



which, taking into account Eqs. (2.3) and (2.5), can be written as 
E 



.dWUt) 
a(t + r) 



dt 

Finally (cf. Eq. (2.4)) 



ft+T 

pk I e- a ^ +T -^5{t-t')dt'. (3.16) 

J —OO 



E [a(t + r)d Wl (t)] = | ^>S' (3.17) 

Therefore, for correlated SV processes future volatility is correlated with past return 
variations, although past volatility and future return variations are completely un- 
corrected. Moreover, p determines the sign of the correlation (3.17) because k is 
positive. In the next section we will show how this property is able to reproduce 
the leverage effect observed in real markets. 

We finally note that, despite the existence of correlations, a(t) and dW\(t) are 
independent random quantities. This is a direct consequence of Ito convention 
for stochastic integrals, because process a(t) is independent of its driving noise 
dWi(t) — £i(t)dt. Hence, as r — > 0~ we have 

E [a(t)dW!(t)} = E [a(t)} E [dW^t)} = 0. (3.18) 

Since dX(t) = a(t)dWi(t). Eq. (3.18) implies that 

E[dX(t)}=0, (3.19) 

in accordance with the name "zero- mean return" given to X{t). 

3.2. Mean reversion 



The effect of mean-reversion refers to the existence of a normal level of volatility 
to which volatility will eventually return. This effect can be observed in financial 
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markets 12 . Practitioners believe that the current volatility is high or low compared 
to a normal level of volatility and they assume that in the long run, forecasts of 
the volatility should all converge to the same normal level. Hence, the average 
of the instantaneous volatility (3.14) should converge to the normal level as time 
tends to infinity, this is done, for instance, by Engle and Patton 5 who define the 
normal level of volatility as the stationary average of the square of instantaneous 
volatility (3.14): 



lim E 

t—fOO 



dX(tf 



dt 



normal level of volatility. 



Note that the limit over time t indicates that process has begun in the infinite past 
and, therefore, the volatility process is in the stationary state. 

Our O-U volatility process appears to be an adequate candidate for describ- 
ing this effect. Let us show this. From Eq. (2.7) and taking in to account the 
independence of <r(t) and dW\(t) and that E [dWf\ = dt, we get 



lim E \dX(t) 2 \a ] = lim E a 2 (t) 



O"0 



dt, 



but the limit t — > oo indicates that volatility has reached the stationary state. 
Hence, the second moment of sigma is given by Eq. (3.10), and therefore 

k 2 

O-U normal level of volatility = 8 2 + — . (3.20) 

Observe that the O-U volatility has a constant and non zero normal level of volatility 
and this is in accordance with the observed mean-reverting property mentioned 
above. 

Moreover, the average over a 2 given by Eq. (3.8), when the volatility is not yet 
in the stationary state, is 

E [a 2 (t)\a ,0] = [a e- at + 9(1 e~ at )] 2 + ^ (l - e^) . 

We observe that this average quickly tends to the normal level (3.20) as at increases. 
This is the reason why the magnitude of a allows us to classify the SV models 
into: (i) fast mean reverting processes when l/a<t, and (ii) slow mean reverting 
processes when l/a>t. 



4. The leverage effect 

It is commonly known that positive or negative sudden changes in the return have 
not the same impact on the volatility. Fisher Black in 1976 13 was the first to find 
empirical evidence on this and observed that the volatility is negatively correlated 
with return variations. A qualitative explanation of this effect is that a fall in the 
stock prices implies an increase of the leverage of companies, which in turn entails 
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more uncertainty and hence higher volatility. Nevertheless, it has also been argued 
that the leverage alone is too small to explain empirical asymmetries in prices 12 . 
Another possible explanation is that news on any increase of volatility reduce the 
demand of stock shares because of investor's risk aversion. The consequent decline 
in stock prices is followed by an increment of the volatility as initially forecasted by 
news, and so on 12 . 

Although the mechanism still lacks of a clear explanation, the leverage effect 
denomination indicates this negative correlation. To our knowledge, this effect has 
only been studied in a qualitative manner until very recently when Bouchaud etal. 
have performed a complete empirical analysis containing new important information 
on this issue. 

Following Bouchaud etal. 11 , we quantify the leverage effect by means of the 
leverage correlation function that we define in the form 

E[dX(t + rfdX(t)] 
C{T) = Var[dX(t)] 2 ' (4 - 21) 

where X(t) is the zero-mean return defined in Eq. (2.6). Bouchaud etal. 11 have 
analyzed a large amount of daily relative changes for either market indices and stock 
share prices and find that 

rM I ~ Ae ~ bT ifT>0 ' uv>\ 

£(t) = \0 ifr<0, (422) 

(A > 0). Hence, there is a negative correlation with an exponential time decay 
between future volatility and past returns changes but no correlation is found be- 
tween past volatility and future price changes. In this way, they provide a sort of 
causality to the leverage effect which, to our knowledge, has never been previously 
mentioned in the literature 10 ' 14 . 

Let us see how our correlated O-U SV model is able to exactly reproduce this 
result. In effect, in the Appendix A we show that 

Vax[dX(t)] = 6> 2 (1 + v 2 )dt (4.23) 

where 

v 2 = k 2 /(2a6 2 ), (4.24) 

is the intensity of volatility fluctuations compared to the expected volatility (see 
Eq. (3.10)). On the other hand we prove in Appendix A that 

f 2pk 6 2 [1 + v 2 er aT ] e- aT dt 2 if t > 0, 
E[dX(t + rfdX(t)] = . ( n (4.25) 

[0 II T < 0. 

Hence the leverage correlation function (4.21) is 



£(r) = 2p 



v\[2a (f + Ve 
(l + ^ 2 ) 2 6> 



for r > 0, (4.26) 
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Fig. 1. The leverage effect in the Dow-Jones index. We plot the leverage function C(t) for the 
Dow-Jones index from 1900 until 2000. We see that there exists a non-negligible correlation when 
t > and negligible when r < 0. Observe that correlation strongly fluctuates when —3 < r < 2. 
We also plot a fit in solid lines with our O-U SV leverage function (4.26). This fit helps us to 
estimate a and p (see Section and Table 2 for more details). 



and 



C(t) =0 for r < 0. 



(4.27) 



We observe that sign of the leverage function is solely determined by the sign of 
p. Therefore, the O-U SV model is able to reproduce the empirically observed 
leverage correlation. In Fig. 1 we show the leverage effect for the Dow- Jones index 
(1900-2000) and plot the leverage function given by Eq. (4.26). 



5. Forecast evaluation 

Any acceptable market model is required to "forecast" the dynamics of the market. 
In other words, the model must be able to reproduce the market behavior and have 
an easy and systematic methodology for estimating its parameters. In our model 
these parameters are p, k, 9, and a . Fouquc etal. for the same model only 
estimate two of these parameters (k and 9) from the empirical second and fourth 
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moment of daily normalized stock changes. Unfortunately they cannot give a clear 
estimation of the volatility auto-correlation time 1/a and of the magnitude of p 
although it was already known that p has to be negative in order to reproduce the 
desired skewness. Let us now see that with the help of the leverage effect we can 
completely estimate all parameters of the model. 

We will perform the empirical estimation of the Dow-Jones index daily return 
changes by approximating dX by AX, i.e., 

dX(t) ~X{t + At) -X{t), 

where At = 1 day. From Eqs. (4.23) and (A. 7) of Appendix A, we have 

Var[AA] = 2 (1 + u 2 )At, Var[AA 2 ] = 2d 4 [4(1 + v 2 ) 2 - 3] At 2 , 

where v 2 = k 2 /(2a0 2 ). From this we get 

1 4 l Var[AX 2 ] 

(1 + ^ 2 ) 2 3 6Var[AA] 2 - [ ' 

Hence, we can estimate the value of v 2 once we know the empirical values of these 
variances. The Dow- Jones index proportionates the daily variances of AX and AX 2 
and, subsequently, gives an estimated value for v 2 . Afterwards, 9 is estimated with 
the knowledge of v 2 and the empirical Var[AX]. In Table 1, we briefly report these 
operations and give the corresponding estimation of v 2 and 6 2 for the Dow- Jones 
index time-series from 1900 until 2000. 

Since (1 + v 2 ) 2 is always positive and (1 -I- j^ 2 ) 2 > 1, we see from Eq. (5.28) that 

VarjAA 2 ] 
S Var[AX]2 < 8 ' 

and the kurtosis, 72 = Var[AA 2 ]/Var[AA] 2 - 2, has the bounds 

< 72 < 6, (5.29) 

showing that the model is never platykurtic. For the Dow- Jones index 72 = 1.72 
and which is consistent with requirement (5.29). However, as Cont reports 16 , there 
exists other markets or even intraday tick data possessing a higher kurtosis outside 
inequality (5.29). 

The leverage effect provides the way for estimating the correlation coefficient p 
and the characteristic time 1/a. Indeed, the best fit of the leverage function (4.26) 
to the Dow- Jones daily data (Fig. 1) gives a characteristic time decay 1/a ~ 20 
days. We thus compare the empirical leverage with our theoretical leverage (4.26) 
when t — > 0^ , 
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Table 1. The O-U SV estimation from the return variances. We estimate the parameters of the 
model from Dow-Jones historical daily returns from 1900 to 2000. We take the variances given by 
Eq. (4.23) and use the identity (5.28) for deriving the estimated quantities v 2 and 9. 



Estimators Dow-Jones daily return Theoretical values 

Var[AX(t)] 1.68 x 10~ 4 9 2 (1 + v 2 )At 

Var [AX(tf] 10.5 x lO" 8 26» 4 [4(^ 2 + l) 2 - 3] At 2 

Parameter estimation V 2 = 0.18 

6 = 18.9% year" 1 / 2 



Table 2. The O-U SV estimation from the leverage. We estimate the parameters p, a and k from 
the fit of the leverage correlation derived from the Dow- Jones stock index data plotted in Fig. 1. 
For doing this, we take the v 2 estimation given by Table 1 and assume that the leverage function 
is given by Eq. (5.30). Observe that magnitudes £(0 + ) and a estimated from the Dow-Jones 
index are of the same order as those given by Bouchaud etal. 11 for a combination of several stock 
indices. 



Estimators Dow-Jones data estimation 

£(0+) -12.5 
a 0.05 day -1 

1/a 19.6 days 

p -0.58 



k = V2ais 2 6 2 1.4 x lO^'Mays" 1 
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and get an estimation of p once we know a = 0.04 days -1 and v 2 = 0.18. Finally, 
we can derive k from the definition given by Eq. (4.24) (sec Table 1). All these 
operations are summarized in Tabic 2. 

In Fig. 2 we simulate the O-U SV resulting process with the parameters es- 
timated above. We follow the random dynamics for AX(t) and compare it with 
the empirical Dow- Jones time series during approximately one trading year. We 
also simulate the ordinary geometric Brownian motion, i.e., assuming a constant 
volatility whose value is directly estimated from data as in Table 1. We clearly 
see in Fig. 2 that the ordinary geometric Brownian model cannot describe neither 
largest nor smallest fluctuations of daily returns. Hence, in comparison with the 
ordinary Brownian modelisation, we conclude that the SV model describes a very 
similar trajectory than that of the Dow-Joncs. This is quite remarkable, because 
we have simulated last year's trajectory using all past data of the Dow- Jones index 
with almost equal results to the actual case. This, in turn, shows the stability of 
parameters. We may thus say that the model is fairly useful for "predicting" the 
stock dynamics (at least for one year period) using market history. 

6. The probability distribution 

We will now obtain the probability distribution of the model. This problem has been 
recently addressed by Shobel and Zhu 17 who, using the Feynmann-Kac functional, 
end up with an expression for the two-dimensional characteristic function of the 
joint process (R(t),a(t)). Here we take a different path that, besides being simpler, 
allows us to get an analytical expression of the return characteristic function, which 
has more practical interest than the joint density. The analysis can be done in terms 
of the return R(t) = In S(t)/ Sq but we prefer to deal with the zero- mean return 
X(t) since, although the calculation is basically identical, the expressions derived 
are shorter and handier. 



6.1. The characteristic function 

Let p 2 {x, a, t\x , (To, h) be the joint probability density of the two-dimensional diffu- 
sion process (X(t), a(t)) described by the pair of SDE's given by Eqs. (2.7) and (2.2). 
This density obeys the following backward Fokker-Planck equation 18 : 



dp 2 , a ,dp 2 1 2 d 2 p 2 , d 2 p 2 l, 2 d 2 P2 
ot octq 2 OTq octoOxq 2 octq 



<*(<to - 6)^ - ^al— - P^o^r^- - 77 fc2 !7Zr> ■ 
with final condition 

p 2 {x, cr, t\x , <t , i) = 5(x - x ) S(<7 - (t )- (6.32) 

Before proceeding further we note that, since all coefficients in Eq. (6.31) are in- 
dependent of x,x ,t and to and the final condition only depends on the differences 
x — x and t — to, then (X(t), <r(t)) is an homogenous process in time and return. 
Therefore, 

p 2 {x,a,t\x ,a ,t ) = p 2 (x - x ,a,t - t \a ). (6.33) 
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Fig. 2. Path simulation and Dow-Jones historical time-series. We show a Dow-Jones daily returns 
sample path (top figure), the O-U SV process simulation, and the geometric Brownian process 
(with constant volatility) simulation (bottom figure). For the OU-SV process, we have taken 
the parameters given by Tables 1 and 2. For the ordinary geometric Brownian motion, a 2 = 
Var[AX]/At and this variance is given in Table 1. The dynamics is traced over approximately a 
trading year (the empirical path approximately corresponds to 1999 trading year). 
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Moreover, one easily sees that the marginal density of the return, 



pxix — xo,t- t \a ) 



p 2 (x - xo,cr,t - t \ao)da, 



also obeys the same partial differential equation than p 2 , Eq. (6.31), which due to 
homogeneity can be written in the form* 



dp x 



= -a(a - 9) 



dp x 



1 o&px 



- pka 



d 2 p x 



da^dx 2 da 2 



(6.34) 



dt ~ v ~ u dn Q ' 2 dx 2 

where px — Px(x — x , t — to\ao) and the initial condition is 

p x (x - x o ,0\a o ) = S(x - x ). 

Partial differential equation (6.34) is the starting point of our analysis. Again, due 
to homogeneity we can assume, without loss of generality that x — and t = 0. 
Observe that Eq. (6.34) is quite involved because of the correlation between the 
volatility and the return, i.e., because of the crossed derivative term. Therefore, 
there seems to be a tremendous task if one tries to attack the problem directly 
from Eq. (6.34). However, Fourier analysis proportionates the necessary tools to 
obtain an analytic solution of the problem. Solution that is expressed in terms of 
the return characteristic function (cf) defined by 



ip x {uJ,t\a ) 



-f 



e lu ' x px(x 1 t\a )dx. 



This is done in Appendix B where we prove that 

tpx(u,t\a ) = exp[- A(u, t)al - B(ui,t)a - C(u,t)], 



where 



A(oj,t) 
B(u,t) 



sinh rjt 



2 y -q cosh rjt + ( sinh rjt 
co 2 a9 ( coshrjt — 1 



rj \ rj cosh rjt + ( sinh rjt 



(6.35) 

(6.36) 
(6.37) 

(6.38) 
(6.39) 

*Note that Eq. (6.31) is the backward Fokkcr-Planck equation while Eq. (6.34) is a forward equation 
for x and t but not for a. 



+ iujpk — a 
(Lua6) 



If ( 
1 12 H — In cosh rjt ^ — sinh rjt 
2 V V 

2£(cosh rjt — 1) + rysinhfyi 



2rji 



rj cosh rjt + Q sinh rjt 



and 



rj = \/ a 2 — 2ipkatu + (1 — p 2 )k 2 lj 2 , Q = a — iujpk. 
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Furthermore, we obtain the unconditional characteristic function, ipx(u>,t), if we 
average over ctq which we assume is in the stationary regime. We thus write 



<px{u,t) = J tpx(u,t\<7 )p st (a )da 
with the stationary pdf given by Eq. (3.12): 

a(a - 6) 



Pst(cr) = 



sj nk 2 /a 



exp 



k 2 



Then from Eqs. (6.35) and (6.40)-(6.41), we get 



y/1 + k 2 A/a 



exp 



-C + 



B 2 k 2 /a - A9B - A6 2 A 
4(1 + k 2 A/a) 



(6.40) 



(6.41) 



(6.42) 



Note that this solution has the right limit when volatility is constant and non- 
random. Indeed, in such a case k = and from Eqs. (6.36)-(6.38) and (6.42) we 
have 



tpx(u,t) = e 



(6.43) 



which is the cf of the zero- mean return when X(t) follows a one-dimensional dif- 
fusive process with constant volatility a = 9. Hence, solution (6.42) appears to be 
consistent with the geometric Brownian motion model. 

6.2. Convergence to the Gaussian distribution 

Let us see that, as t — > oo, the marginal distribution (char2) approaches to the 
Gaussian density. In other words, we will prove a Central Limit Theorem for the 
model. The starting point of this analysis is Eq. (6.42) that as at > 1 can be 
written in the following simpler form 



yx{u,t) ~ exp 



rf 



iujpk 



V)t/2 



(at > 1). 



(6.44) 



This is not a Gaussian distribution yet, since 77 defined by Eq. (6.39) is an irrational 
function of lo. We have to assume an extra requirement. Specifically, we suppose 
that 

- < 1, (6.45) 

a 

which means that volatility is weakly random. Indeed, k is the strength of the 
volatility driving noise (the so-called volatility of volatility) while a tell us how 
large is its deterministic drift (see Eq. (2.2)). Therefore the ratio k/a measures in 
some way the degree of volatility randomness. Taking into account Eq. (6.45) we 
write (cf. Eq. (6.39)) 



i] = a 



k 
)— 

a 



a 2 



iUp— + TT^ 2 ^ + O I 
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Substituting this into Eq. (6.44) yields 

V >x(u,t)~cxp{-uj 2 [l + is 2 + 0(k/a)]6 2 t/2}, (at > 1), (6.46) 

where v 2 = k 2 /2a8 2 . The Gaussian density (6.46) proves the Central Limit Theo- 
rem in our case. 

6.3. Cumulants 

Cumulants are defined as follows 

d n 

Kn = a^fo* (<«>.*)] 

and are very useful for deriving statistical properties of the model. For instance, 
the second cumulant reads 

k 2 = 2 (l + v 2 ) t, (6.47) 

which results to be the integrated variance of dX(t) (cf. Eq. (A. 7)). The third and 
fourth cumulants are respectively 



(6.48) 



+p 2 (1 - e- 2at ) 



.3 = {2 [at (1 - e-*)] + y [2a* - (l - e"**)] } , 

and 

2at + atp 2 (6 + 4e~ at ) - (2 + 12p 2 ) (l - e" at ) 

v 2 [2at + 8atp 2 (l + e" 2 " 4 ) - (l - e" 2 "*) (l + 8p 2 )] j. (6.49) 

In terms of cumulants kurtosis is given by 

72 = 5> (6-50) 

and measures the tails of the distribution compared to the Gaussian distribution. 
The kurtosis can be thus obtained using the second and fourth cumulants given by 
Eqs. (6.47) and (6.49). The resulting expression is very similar to that of the fourth 
cumulant (6.49) but with an extra constant factor. Its asymptotic limits are rather 
simple to derive and read 

<-«»■ <™> 

and 

6^> 2 (l + 4p 2 )+4(l + p 2 )] 1 
72 7^)2 ~ t («* » !)• (6-52) 
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From Eq. (6.51), we observe that even dealing with an infinitesimal time there will 
exist a non negligible kurtosis. Conversely, from Eq. (6.52) we see that kurtosis goes 
to zero as time increases and that the convergence is slow going as 1/t. In addition, 
we observe that for short times kurtosis does not contain the correlation coefficient 
p but in the long run a non zero p magnifies the kurtosis of the distribution (c/. 
Eqs. (6.51)-(6.52)). 

Let us now turn our attention to skewness 71 defined in terms of third and 
second cumulant as 

7i = -^75 (6-53) 

and that quantifies the byass in the return distribution. A negative skewness in- 
dicates that returns are more likely to decrease than to increase, while a positive 
skewness indicates a higher probability of a return raising than a decline. Empirical 
observations have found that financial markets have an slightly negative skewness 16 . 

Similarly to the kurtosis derivation, the skewness defined by Eq. (6.53) is ob- 
tained with Eqs. (6.47)-(6.48) and the resulting expression is very similar to that of 
Eq. (6.48) but with an extra constant factor. We limit ourselves to the asymptotic 
cases 

71 - 3p V V2al (for at < 1), (6.54) 
V^ 2 + 1 



and 



7i ~ 6 P Pr -7= ( for at > !)• ( 6 - 55 ) 



( V 2 + lf 2 V2al 

From these equations we see that both at short and long times, skewness vanishes. 
However, it decreases very slowly, more slowly than kurtosis (compare Eqs. (6.52) 
and (6.55)). Finally, we note that skewness is proportional to p and, in consequence, 
the sign of p not only determines the leverage correlation sign but also the skewness 
sign. Empirical observations of leverage and skewness indicate that p must be 
negative 16 . 



6.4. Tails 

It is well established that distribution of prices have heavy tails. There exists several 
empirical studies quantifying this fact (see for instance Mantegna and Stanley 19 or 
Plerou etal. 20 . Let us now study the existence of fat tails in our SV model. Recall 
first that for long times, at^>l, the probability distribution is practically Gaussian 
and there is no fat tail to look for. Therefore, we will search for heavy tails at small 
to moderate times, i.e., when Central Limit Theorem is not applicable. 

The tails of the distribution are determined by the shape of the density function 
Px{x, t) when x is large. A well known fact from Fourier analysis is that the large 
x behavior of px{x, t) is given by the small lj behavior of its characteristic function 
<px(u>,t) 21 . Therefore, tails are derived from the characteristic function (6.42) by 
only keeping the first orders in u. 



17 



When u> is small but time is not too long (i.e., at <~ 1) the expressions for 
B(u>,t) and C(ui,t) are approximately given by (c/. Eqs. (6.36)-(6.38)) 



A(o;,t)~-(l-e- 2at ), 



and 



C(w, t) - (u 2 8 2 + iujpk - a)t/2 - (6» 2 /4a)[2(l 



l ) 2 + l-e- 2at ]u 2 



+- ln[l - ipfciw + (fc74a 2 )(2ai - 2a 2 t A p A - 1 + e~ 2at )uj 2 }. 
Thus from Eq. (6.42) we have 

[1 + (fc 2 /4a 2 )(l - e ~2at)]-i/2 exp (-^ 2 g 2 t/2 - ia^/2) 
Px(w,t) ~ _ + ( A .2/4 a 2)( 2a t _ 2a 2 t 2 p 2 - 1 + e- 2 « t )w 2 ] 1 /2 ' 

Again, taking into account that u> is small and t is moderate we get 

1 



1 - 2ia(t)w + b(t)uj 2 ' 



where 



a(i) = p/rf/4, and &(i) = fc 2 i(2 - p 2 at)/8a 
The inverse Fourier transform for this asymptotic cf is 

~ [^a(t) 2 +b(t)\x\ - a{t)x) 



(6.56) 
(6.57) 



Px(x,t) 



exp 



\x\ — > oo). 



Hence the tails of the zero-mean return have an asymmetric exponential decay given 

by 



Px(x,t) 



1 



exp 



L (V^) 2 + W) - aif)) x 



b(t) 



(x — » oo) 



(6.58) 



and 



Px(x,t) 



1 

Va(i) 2 +6(t) 



exp 



^ ( V«W 2 + &(*) + «(*)) s 



(.x — » — oo). 



(6.59) 

Since a(i) = pkt/A, we see that the sign of p will determine which is the heaviest 
tail. When p is negative the fattest tail is the one representing losses and when 
p > the fattest tail corresponds to profits. If p = there is no difference between 
the two tails. 

Finally, let us guess how the tails of price (not return) distribution are. We first 
recall that the asymptotic expressions (6.58)-(6.59) refer to the marginal distribution 
of X(t), that is, regardless the value of volatility at time t and after averaging over 
the initial volatility. In order to obtain the asymptotic form of price distribution 
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Ps(S,t) out of the asymptotic form of px(x,t) we must know what is the relation 
between S(t) and X(t). For the general case this relation is given by Eq. (A.l) of 
Appendix A and, since it corresponds to the two-dimensional case when no average 
and marginal distribution have been performed, Eq. (A.l) involves X(t), S(t) and 
cr(t). We conjecture that if time is not too large X(t) ~ \n[S(t)/ So]. Therefore, for 
the price distribution ps{S,t) we have the following power laws: 



Ps(S,t) 



' (S-0), and Ps (S,t)~— Ij-r (S^oo) (6.60) 



where 



y±{t) = 1 + ^ [±^a(ty+b(t)-a(t) 



7. Conclusions 

The stochastic volatility (SV) models are a possible way out to the observed in- 
consistencies between the geometric Brownian model and real markets. The SV 
models, as their name indicates, assume the original log-Brownian model but with 
the volatility a being random. From 1987 on, there have appeared several works 
extending and refining stochastic volatility models but most of them are basically 
designed to reproduce empirical option prices. We have assumed that the volatility 
follows one of the simplest stochastic volatility models still showing mean reversion, 
i.e., the Ornstein-Uhlenbeck process, but also allowing for correlations between the 
random fluctuations of the price and volatility processes. 

Recent efforts in the study of the empirical statistical properties of the spec- 
ulative markets has led to a list of stylized facts that all proposed model should 
accomplish. We have shown that the model is able to reproduce these facts. More 
specifically, the volatility process is mean reverting and have found that it is able 
to quantitatively reproduce the recently observed leverage effect. We have also es- 
timated all parameters of the model and observed that a simulated path (using the 
estimated parameters) is very similar to the sample path of the historical evolution 
of the Dow- Jones daily index (1900-2000). Finally, we have derived the character- 
istic function of the process, obtained its kurtosis and skewness, and observed the 
power-law decay for the tails of the price distribution. The results herein derived 
show that the SV models are good candidates for describing not only option prices 
but market dynamics as well. 

We finally point out two further extensions of our work. First, the resulting 
characteristic function with all parameters already estimated can be used to derive 
option prices. Second, the observed volatility is not definitively attached to any 
specific model. There exists a large class of volatility processes and, since we now 
have a systematic way of estimating parameters of the process, a deeper analysis 
on which model gives a better explanation of the stylized facts can also be very 
interesting. Other existent volatility models require a more involved mathematical 
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analysis and it would be also interesting how far we can go using the methods 
presented in this paper. These two points are under present investigation. 
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Appendix A. The zero-mean return 

The zero- mean return has been defined through its differential dX by Eq. (2.6). Let 
us first prove that X{t) is explicitly given by 

X(t) = HS(t)/S Q ] -fa -I J* a 2 (t')dt'. (A.l) 

Indeed, if we apply the Ito lemma to Eq. (A.l) and, as usual, keep orders smaller 
than dt 3 / 2 we have 

dX(t) = dS/S + ^(dS/Sf - [idt - ^a 2 dt, 

but dS/S = [idt + <rdW 1 and {dS/S) 2 = a 2 dt, then we obtain Eq. (2.7): 

dX = adW u 

and this proves the validity of Eq. (A.l). 

We will now derive several averages concerning dX. We know from Eq. (3.19) 
that 

E[dX(t)]=0. (A.2) 
Again taking into account the independence of a(t) and dW\{t) we write 

E [dX 2 ] = E [a 2 ] E [dW^tf] , (A.3) 

but E [dWf] = dt, and using Eq. (3.10) we have 

E [dX 2 ] = (0 2 + k 2 /2a) dt. (A.4) 

As to the fourth moment, 

E [dX 4 ] = E [a 4 ] E [dW 4 ] , 

we note that E [dWf\ = 3E [dW?] 2 = Ut 2 and we evaluate E [a 4 ] using the 
stationary pdf (6.41). Hence 

E [dX 4 ] = 3 [3(fc 2 /2a) 2 + 6(k 2 /2a)0 2 + 4 ] dt 2 . (A.5) 
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The variances of dX 2 and dX are obtained from Eqs. (A. 4) and (A. 5) and read 

Var[dX(t)] = (9 2 + k 2 /2a) dt, (A.6) 
Var[dX(t) 2 } = 2\0 4 + A(k 2 /2a) 2 + 8(k 2 /2a)d 2 ]dt 2 . (A.7) 

We finally derive the following correlation function: 

E[dX(t + T) 2 dX(t)] = E [a{t)dWi{t)o{t + r) 2 dW 1 {t + t) 2 ] . 

Note that all variables are Gaussian which allows us to decompose the rhs of this 
equation into average pairs, taking also into account that dW\{t + r) 2 = dt we can 
write 

E [dX(t + T) 2 dX(t)] = ^2E[j(t)a(t + T)]E[j(t + T)dWi(t)] 

+E [a(t + t) 2 ] E[a{t)dWi{t)]}dt 

Combining this with Eqs. (3.17) and (3.18), we get 

e [dx {t + rfdxit)] = { 2 * E wm* + ^ e ~ aTdt2 

and the volatility correlation (3.11) allows us to write 

f 2pk \9 2 + (k 2 /2a)e- aT ] e - aT dt 2 if r > 0, 
E[dX(t + T) 2 dX(t)}=\ n . t n (A.8) 

if r < 0. 



Appendix B. The characteristic function of return 

In this Appendix we will obtain the expression given by Eq. (6.35) for the marginal 
characteristic function y>x{k>,t\(Jo) of the two-dimensional diffusion process (X,cr) 
whose joint density P2{%, cr, t\xo, <jq, to) is the solution to the final problem posed by 
Eqs. (6.31)-(6.32). 

The marginal characteristic function (cf) of process X(t), 

/CO 
e iojx p x {x,t\a )dx, 
-oo 

allows us to write Eq. (6.34) in the following simpler form: 

depx 1 2 d 2 tp x . r . , , aw d V x 1 2 2 , D1 ^ 

~df = 2 ~d^f [lUJpka ° ~ a{aa ~ 6)] ^- 2°°" * x ' (R1) 

The initial condition is 

tp x (u,0\a ) = l. (B.2) 
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By direct inspection one can easily see that the solution to problem (B.1)-(B.2) 



is 



ipx(tJ,t\a ) = exp[-A(w,t)al - B{uj,t)a - C{u,t)] (B.3) 

where functions A(w,t), B(w,t) and C(u,t) are the solution of the following set of 
ordinary differential equations 

A = -2k 2 A 2 -2{a-iujpk)A + uj 2 /2 (B.4) 
B = — [2k 2 A +(a- iujpk)] B + 2aBA (B.5) 
C = k 2 (A- B 2 /2) + a6B, (B.6) 

with initial conditions 

A(u,0) =B(lj,0) = C(w,0) = 0. (B.7) 

Note that Eq. (B.5) is a linear equation and that the rhs of Eq. (B.6) does not 
contain C(t). Therefore, their formal solutions are straightforward and read 

B(t) = 2a6 f A(t') exp -{a - iupk){t - t') - 2k 2 f A(t")dt" dt' , (B.8) 
Jo I Jt' 

C(t) = k 2 [ [A(t') - B 2 (t')/2] dt' + a6 f B(t')dt'. (B.9) 
Jo Jo 

On the other hand Eq. (B.4) is more involved since it is a Ricatti equation. However, 

the definition of a new dependent variable 

turns Eq. (B.4) into the following linear second-order equation with constant coef- 
ficients: 

y + 2(a — iujpk)y — k 2 uj 2 y = 0. 
The solution to this equation is 

y{t) = C 1 e x + t + C 2 e x -\ 

where are arbitrary constants and 

A± = a — iupk ± \/Ja — iojpk) 2 + k 2 uo 2 . 

Substituting this into Eq. (B.10) yields 

A+e^ + A-CtyCQe*-* 
[ > 2fc 2 [e^* + (C 2 /Ci)e A -*] ' 

Now the initial condition A(0) = gives C2/C1 = — A + /A_ and the substitution of 
A± allows us to write A(t) in the form given by Eq. (6.36). Finally the substitution 
of Eq. (6.36) into Eqs. (B.8)-(B.9) results in Eqs. (6.37)-(6.38). 
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